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ABSTRACT 

Hot subdwarf B stars (sdBs) are evolved, core helium-burning objects located on the extreme horizontal branch. Their formation 
history is still puzzling as the sdB progenitors must lose nearly all of their hydrogen envelope during the red-giant phase. About half 
of the known sdBs are in close binaries with periods from 1.2 h to a few days, a fact that implies they experienced a common-envelope 
phase. Eclipsing hot subdwarf binaries (also called HW Virginis systems) are rare but important objects for determining fundamental 
stellar parameters. Even more significant and uncommon are those binaries containing a pulsating sdB, as the mass can be determined 
independently by asteroseismology. Here we present a first analysis of the eclipsing hot subdwarf binary V2008-1753. The light curve 
shows a total eclipse, a prominent reflection effect, and low-amplitude pulsations with periods from 150 to 180 s. An analysis of the 
light- and radial velocity (RV) curves indicates a mass ratio close to q = 0.146, an RV semi-amplitude of K = 54.6 kms^*, and an 
inclination of i = 86.8°. Combining these results with our spectroscopic determination of the surface gravity, log g = 5.83, the best¬ 
fitting model yields an sdB mass of 0.47Mo and a companion mass of 69Mjup. As the latter mass is below the hydrogen-burning limit, 
V2008-1753 represents the first HW Vir system known consisting of a pulsating sdB and a brown dwarf companion. Consequently, it 
holds great potential for better constraining models of sdB binary evolution and asteroseismology. 

Key words, stars: subdwarfs, binaries: eclipsing, binaries: spectroscopic, stars: brown dwarfs, stars: fundamental parameters, stars: 
individual: V2008-1753, stars: oscillating 


1. Introduction 

Hot subdwarf B (sdB) stars are core-helium burning stars with 
very thin hydrogen envelopes that are found on the extreme hor¬ 
izontal branch (EHB) (see Heber 2009, for a review). While the 
Future evolution of sdB stars is quite certain (they will become 
white dwarfs), their prior evolutionary paths remain to be fully 
resolved. Nonetheless, one thing seems to be assured: the rel¬ 
atively high fraction of sdBs in binaries (Maxted et al. 2002; 
Napiwotzki et al. 2004) implies that binary interactions play a 
crucial role in their formation. Han et al. (see 2002, 2003, and 
references therein) proposed several binary-related mechanisms 
as formation channels for sdB stars: 


- common envelope ejection leading to short-period binaries 
with periods from 0.1-10 days with either white dwarf or 
low mass main sequence companions. From binary popula¬ 
tion synthesis (Han et al. 2002) a mass distribution that peaks 
sharply around 0.47 M© is expected. This mass is called the 
canonical mass. 

- stable Roche lobe overflow resulting in sdB binaries with 
masses around 0.47 Mq and orbital periods of 10-100 d 

- double helium white dwarf mergers giving rise to single sdB 
stars with a wider distribution of masses. 


Common envelope (CE) evolution is believed to be an impor¬ 
tant process for a large number and wide diversity of binary stars. 
In particular, this formation process is vital to compact binary 


systems, as many of them once must have been orders of magni¬ 
tude wider than their present-day separations. Although binary 
population synthesis (BPS) models are, in general, able to pro¬ 
duce compact binaries through common envelope evolution, this 
important phase in binary evolution is still poorly understood 
(see Ivanova et al. 2013, for a review). Many of the CE-related 
parameters in BPS codes remain unconstrained by observations 
(Clausen et al. 2012), including common-envelope ejection effi¬ 
ciency, minimum mass for core He ignition, and envelope bind¬ 
ing energy. Consequently, the investigation of post-CE systems, 
especially those containing pulsators or exhibiting eclipses, can 
improve our understanding of this highly-important stage of 
stellar evolution. 

Some have proposed that a planet or brown dwarf com¬ 
panion might also be able to trigger the loss of envelope mass 
in the red-giant phase of the sdB progenitor (Soker 1998; 
Nelemans & Tauris 1998; Han et al. 2012). When the host star 
evolves and becomes a red giant, close substellar companions 
will be engulfed in a common envelope. The outcome of this in¬ 
teraction is unclear, but there are three distinct possibilities: the 
substellar object can either survive, merge with the core of the 
red giant, or evaporate. The latter two outcomes might provide 
an explanation for some of the observed single sdBs. Discover¬ 
ies of brown dwarfs in close, eclipsing sdB systems (Geier et al. 
2011; Schaffenroth et al. 2014b) support the idea that substel¬ 
lar companions are sufficient for ejecting the envelopes of red 
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giants. This scenario is supported also by the discovery of a 
brown dwarf (Mbd = 0.053 + 0.006 Mq) in a close orbit (0.08 d) 
around a low-mass white dwarf (Maxted et al. 2006). Conse¬ 
quently, substellar companions are able to influence late stellar 
evolution. 

Eclipsing binaries are of special interest as they allow the de¬ 
termination of the masses and radii of both components as well 
as the period and separation of the system. The mass in particu¬ 
lar is vital to constraining sdB formation models. Eclipsing sdB 
binaries with low mass main sequence or brown dwarf compan¬ 
ions (HW Vir systems) are easily recognizable by the shapes of 
their light curves, which are dominated by both eclipses and a 
strong reflection effect due to the large temperature difference 
between both components. The amplitude of the reflection ef¬ 
fect depends on the separation distance, the temperature of the 
sdB, and the albedo of the companion. In addition to shedding 
light on sdB evolution, studies of these relatively rare, post¬ 
common envelope systems can be used to constrain current mod¬ 
els of common-envelope evolution (e.g. Zorotovic et al. 2010, 
for eclipsing WD binaries). 

Additional information concerning the structure and evolu¬ 
tion of sdBs may be found by studying their pulsations. The 
sdBVi stars, discovered by Kilkenny et al. (1997) and indepen¬ 
dently theoretically predicted by Charpinet et al. (1996), are low 
amplitude, multimode pulsators with typical periods ranging be¬ 
tween 80-600 s. Their pulsation amplitudes are generally of the 
order of a few mmag. The short periods, being of the order of and 
shorter than the radial fundamental mode for these stars, sug¬ 
gest that the observed modes are low-order, low-degree p-modes 
(Charpinet et al. 2000). The known sdBVr stars occupy a region 
in the Teff - log.g plane with effective temperatures between 
28 000 K and 36 000 K and surface gravities (log g) between 
5.2 and 6.2. Only 10% of all stars falling in this region show 
pulsations. Moreover, also sdBs with slow high radial-order g- 
mode oscillations with periods on the order of 2000 - 8000 s 
and termperatures of 22000 - 29000 K were found (Green et al. 
2003). In the overlap region also hybrid pulsators showing p- 
mode as well as g-mode pulsations exist (Schuh et al. 2006). 
Asteroseismology of these stars allows for an independent and 
accurate determination of the sdB mass. 

The problem with typical HW Vir systems has been 
that they are single-lined spectroscopic binaries (see e.g. 
Schaffenroth et al. 2013, 2014b), and as such it is normally not 
possible to derive a unique mass ratio. HW Vir systems har¬ 
boring a pulsating sdB primary, however, offer additional pos¬ 
sibilities, as the stellar properties can be constrained by the light 
curve and asteroseismological analyses. Until recently, only two 
HW Vir systems with pulsating sdBs were known. The first such 
object - NY Vir - was found to be a pulsating sdB in an eclipsing 
binary by Kilkenny et al. (1998). It shows more than 20 pulsation 
modes with amplitudes of several mmag. An astroseismic analy¬ 
sis was able to determine the stellar parameters of this system 
independently from the lightcurve analysis (Van Grootel et al. 
2013). 0stensen et al. (2010) found another pulsating sdB -)-dM 
HW Vir binary (2M1938H-4603). Unfortunately, the amplitudes 
of the pulsations, which were detected by Kepler are so small 
that they cannot be observed by ground-based telescopes. Thus, 
this star is therefore not an ideal target for asteroseismological 
modelling. 

Van Noord et al. (2013) reported the discovery of a promis¬ 
ing new HW Virginis system, V2008-1753 (CV=16.8 mag), 
which was found during an automatic search for variable stars 
conducted with the 0.4 m Calvin College Robotic Telescope 
in Rehoboth, New Mexico. Their relatively noisy light curve 


showed eclipses and a strong reflection effect. Interestingly, this 
sdB binary has an orbital period of only 1.58 h, the shortest pe¬ 
riod ever found in a HW Virginis system. Here we present the 
first thorough analysis of this unique system, along with the dis¬ 
covery of low-amplitude pulsations in the sdB primary. Section 
2 describes the observational data. The analysis is explained in 
Sects. 3 (spectroscopy) and 4 (photometry). Evidence for the 
brown dwarf nature of the companion is discussed in Sect. 5. 
Einally, we end up with conclusions and suggest further oppor¬ 
tunities that are offered by this one-of-a-kind system. 

2. Observations 

2.1. Time-Series Spectroscopy 

We used the Goodman Spectrograph on the 4.1-m SOuthern 
Astrophysical Research (SOAR) telescope to obtain time-series 
spectroscopy of V2008-1753 over a full orbital cycle and deter¬ 
mine the orbital velocity of its primary sdB component. A 1.35” 
longslit and a 930 mm ' VPH grating from Syzygy Optics, EEC, 
were employed to cover the spectral range 3600-5250 A with an 
approximate resolution of 3.8 A (0.84 A per binned pixel). The 
position angle was set to 197.8 degrees E of N so that we could 
place a nearby comparison star on the slit, with the intention of 
characterizing and removing instrumental flexure effects. In or¬ 
der to maximize our duty cycle, we binned the spectral images 
by two in both the spatial and dispersion directions and read out 
only a 2071 x 550 (binned pixels) subsection of the chip. Each 
exposure had an integration time of 90 s, with 8 s of overhead 
between successive images, yielding a duty cycle of approxi¬ 
mately 92%. Overall, we acquired 70 spectra of V2008-1753 be¬ 
tween 23:49:20.87 UT (2013-08-31) and 01:42:07.17 UT (2013- 
09-01). The airmass decreased from 1.22 to 1.03 over this time 
period. Upon the completion of the time series, we obtained sev¬ 
eral EeAr comparison lamp spectra and quartz-lamp spectra for 
wavelength calibration and flat-fielding. 

Standard routines in IRAE, primarily ccdproc, were used to 
bias-substract and flat-held the spectral images. Given the low 
thermal noise in the spectrograph system, we did not subtract 
any dark frames as we wanted to avoid adding noise to the im¬ 
ages. We used apall to optimally extract one-dimensional spec¬ 
tra and subtract a fit to the sky background for both the target 
and constant comparison star. Einally, we wavelength-calibrated 
each spectrum using the master EeAr comparison lamp spectrum 
taken at the end of the time-series run. The resulting individual 
spectra of V2008-1753 had a signal-to-noise ratio (S/N) around 
15-20 pixel ', while the individual spectra of the comparison 
star (which looked to be G or K-type), had a S/N twice as high. 

2.2. Time-Series Photometry 

High-precision photometry was acquired with SOAR/Goodman 
through i’ and g’ filters on 15/16 August 2013. Although our 
primary goal was to model the binary light curve, our secondary 
goal was to look for stellar pulsations, and thus we used an in¬ 
strumental setup that had both a high duty cycle and a Nyquist 
frequency above those of most known sdB pulsation modes. We 
binned the images 2x2 and read out only a small 410 x 540 
(binned pixel) subsection of the chip that included both the target 
and 10 comparisons stars with the same approximate brightness 
level. We achieved an 87% duty cycle with 25-s exposures for the 
g’ light curve, and a 81.4% duty cycle with 14-s exposures for 
the i’ light curve. In both cases, we observed the field for at least 
one full orbital cycle near an airmass of 1.1. We concluded each 
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night with a set of bias frames and dome flats. Thermal count 
rates were too low to warrant the acquisition of dark frames. 

All object images were bias-subtracted and flat-fielded us¬ 
ing IRAF’s ccdproc package. We extracted light curves with an 
IDL program we wrote that uses the function APER, which is 
based on DAOPHOT (Stetson 1987). We produced light curves 
over a wide range of aperture radii and selected the aperture that 
maximized the S/N in the light curve. To mitigate the effects of 
atmospheric extinction and transparency variations, we divided 
the light curve of V2008-1753 by the average of those of the 
constant comparison stars, after verifying they were indeed non¬ 
variable. Residual extinction effects are often removed in light 
curves of this duration by fitting and dividing the light curve by 
the best-fitting parabola. However, given the large-amplitude bi¬ 
nary signals present, this process would distort the actual stellar 
signal. Instead, we fit straight lines through points of the same 
phase and, informed by these fits, removed the overall ’tilt’ in 
each light curve. While not perfect, this process helps to mollify 
the effects of residual extinction variations. We then divided each 
light curve by its mean value and subtracted a value of one from 
all points to put them in terms of fractional amplitude variations. 

3. Spectroscopic analysis 

3.1. Radial velocity curve 

Radial velocity shifts were determined by measuring the posi¬ 
tions of the hydrogen Balmer profiles H/? through H9; although 
HIO, Hll, and several He I lines were also present, they were 
too noisy to provide reliable positions. We used the MPFIT rou¬ 
tine in IDL (Markwardt 2009), which relies on the Levenberg- 
Marquardt method, to fit simple inverse Gaussians to the line 
profiles and determine their centroids. The only available guide 
star near our field was significantly redder than the sdB target, 
and, consequently, its use led to a gradual shift in the slit align¬ 
ment over the course of our observations (Goodman had no at¬ 
mospheric dispersion corrector at the time); this drifting results 
in a colour-dependent velocity shift. Additionally, instrumen¬ 
tal flexure as the Nasmyth cage rotates also affected the stars’ 
alignment on the slit, although only slightly. We removed both 
of these time-dependent wavelength solution effects by track¬ 
ing the absorption-line features in the constant comparison star. 
Figure 1 presents the resulting radial velocity curve for V2008- 
1753, plotted twice for better visualisation. 

We again used MPFIT to fit a sine wave to the data and de¬ 
termine the semi-amplitude of the velocity variation; the orbital 
period and phase were fixed during this process. We derive an 
orbital velocity of /T = 54.6 + 2.4 km s ' for the sdB primary. 
Eccentric orbits were also fitted to the radial velocity curve, but 
as we currently have no reason to prefer them over e = 0, we 
continue the analysis under the assumption that the orbit is cir¬ 
cular. Residuals from the best-fitting sine wave are shown in the 
bottom panel of Figure 1 and are consistent with noise. The mean 
noise level in the Fourier transform of the residual velocity curve 
is 2kms“*. 

3.2. Atmospheric parameters 

In preparation for determining the atmospheric parameters of the 
sdB, we first de-shifted all individual spectra according to our 
orbital solution above and then co-added them to improve the 
overall S/N. We fit synthetic spectra, which were calculated us¬ 
ing local thermodynamical equilibrium model atmospheres with 
solar metallicity and metal line blanketing (Heber et al. 2000), to 



0 0.25 0.5 0.75 1 1.25 1.5 1.75 2 

orbital phase 


Fig. 1. Top panel: Radial velocity curve for the sdB primary in V2008- 
1753, plotted twice for better visualization. The solid line denotes the 
best-fitting circular orbit to the data. Bottom panel: Residuals after sub¬ 
tracting the best-fitting sine wave from the data. 

the Balmer and helium lines of the co-added SOAR spectrum us¬ 
ing SPAS (SPectral Analysis Software, Hirsch 2009). The best¬ 
fitting synthetic spectrum had the following orbital parameters; 

Teff = 32800 + 250 
logg = 5.83 + 0.04 
logy = -2.27 + 0.13 

with 1-cr statistical errors determined by bootstrapping. For 
some sdB systems with reflection effects, an analysis of spectra 
with sufficiently high signal-to-noise taken at different phases 
shows that the atmospheric parameters apparently vary with 
phase (e.g. Schaffenroth et al. 2013, 2014b). These variations 
can be explained by the companion’s contribution to the spec¬ 
trum (reflection only) varying with orbital phase. Systems with 
similar parameters show such variations in temperature and sur¬ 
face gravity on the order of 1000 - 1500 K and 0.1 dex, re¬ 
spectively. To account for the apparent change in the parame¬ 
ters we formally adopt the values determined from the co-added 
spectrum, which represents a mean value, with a larger error; 
Teff = 32800 ± 750 K and log g = 5.83+0.05 for the sdB. Figure 
2 shows the corresponding fit of the Balmer and helium lines. 
We excluded He from the fit, as this line is mostly blended with 
the Ca II H-line and hence is less well represented by this fit. 

The Teff - log g diagram is displayed in Fig. 3 and shows that 
V2008-1753 lies in the middle of the extreme horizontal branch. 
Although it was previously suggested that HW Virginis systems 
cluster together only in a small part of the - log g diagram 
(Schaffenroth et al. 2014b), the position of V2008-1753 seems 
to go against this hypothesis. However, it is still apparent that 
most of the known HW Vir systems and reflection effect binaries 
(both sdB+dM systems with different inclinations) concentrate 
in a distinct region between a of 26000-30000 K and a log g 
of about 5.3 to 5.7, only at the edge of the instability strip. There 
are five exceptions of binaries with sdOB star primaries and M 
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Fig. 2. Fit of the Balmer and helium lines in the co-added SOAR spec¬ 
trum. The solid line shows the measurement, and the dashed line shows 
the best fitting synthetic spectrum. As the H e seems to be effected by a 
blend of the Ca line next to it, it was excluded from the fit. 


star companions at higher temperatures, which are possibly just 
more evolved. Yet, the three HW Vir systems showing short- 
period p-mode pulsations with amplitudes observable from the 
ground lie in a different part of the - log g diagram, nearer to 
the He-MS, in the central part of the instability strip for sdBViS, 
as expected. In contrast to that the other sdB binaries with either 
white dwarf secondaries or companions of unknown type do not 
show any clustering but are uniformly distributed over the EHB. 
This could indicate that the sdBs with low-mass main sequence 
companions differ from the sdBs with white dwarf companions. 


4. Photometric analysis 

4.1. Pulsations 

Both the g’ and i’ light curves from the SOAR telescope ex¬ 
hibit pulsations too low in amplitude to have been detected in the 
data analyzed by Van Noord et al. (2013). In order to disentan¬ 
gle the pulsations from the binary effects in the light curve, we 
used an approach similar to that demonstrated by Vuckovic et al. 
(2007a). First, we fit the eclipses and reflection effect as de¬ 
scribed in Section 4.2) and subtracted the best-fitting model from 
the observed light curve. The original g’ and i’ filter light curves 
are displayed in Fig. 4, along with the same curves after the sub¬ 
traction of the binary signal. The pulsations are visible by eye 
in the g’ curve but less apparent in the i’ data, due to its lower 
S/N and lower pulsation amplitudes at redder wavelengths. The 
large-scale trends in both residual light curves are likely due to 
residual atmospheric extinction and transparency variations. 

We calculated the Fourier transformations (FTs) of the light 
curves using FAMIAS' (Zima 2008). The resulting FTs are dis¬ 
played in Fig. 5. We detect at least four independent pulsation 
modes with periods ranging from 2.5 to 3 min and amplitudes < 
4 ppt. The best-fit frequencies, along with their amplitudes and 
S/N, are listed in Table 1. The mode with the second-highest 

* http://v™w.ster.kuleuven.be/~zima/fainias 
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Fig. 3. Teff - logg diagram of the HW Vir systems. The helium 
main sequence (HeMS) and the EHB band (limited by the zero-age 
EHB, ZAEHB, and the terminal-age EHB, TAEHB) are superimposed 
by evolutionary tracks by Dorman et al. (1993) for sdB masses of 
0.471, 0.473, and 0.475 Mq. The positions of the HW Vir systems 
with pulsating sdBs - V2008-1753 (this work, with error bars), NY 
Vir (Van Grootel et al. 2013), 2M1938-H4603 (0stensen et al. 2010), 
and PTFl J072456-I-125301 (Schindewolf et al. submitted, Kupfer 
priv. com.) - are marked by red squares. Blue dots mark the posi¬ 
tions of other HW Vir-like systems (Drechsel et al. 2001; Foretal. 
2010; Maxted et al. 2002; Klepp & Rauch 2011; 0stensen et al. 
2008; Wood&Salfer 1999; Almeida et al. 2012; Barlow et al. 2013; 
Schaffenroth et al. 2013, 2014a, Schaffenroth et al., in prep, Kupfer, 
priv. com.). The positions of the two HW Vir systems with BD com¬ 
panions J1622 (Schaffenroth et al. 2014b) and J0820 (Geier et al. 2011) 
are indicated by the blue dots with error bars. The black triangles 
mark sdB-tdM systems showing a reflection effect but no eclipses 
(Kupfer et al. 2015, and references therein). The green, open dots repre¬ 
sent other sdB-t-WD binaries or sdB binaries with unknown companion 
type (Kupfer et al. 2015). The approximate location of the sdBVj insta¬ 
bility strip is indicated by an ellipse. 



Fig. 4. g' (dashed line) and i’ (solid line) filter light curves of V2008- 
1753 taken with SOAR. The two light curves were taken in subsequent 
orbital cycles. The sub-figure shows the pulsation signal after the sub¬ 
traction of the binary signal by the best-fitting light curve model. 
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frequency [Hz] frequency [Hz] 

Fig. 5. Fourier transform of the pulsation signal shown in the sub-figure of Fig. 4. The left and right figures show the FTs of the g’ and the i’ filter 
light curves, respectively. 


Table 1. Pulsation frequencies and amplitudes 



^a, b 

[mHz] 

amplitude 

[PPt] 

phase"^ 

S/N 

El 

6.565+0.005 

3.5+0.2 

0.739+0.019 

14.5 

E2 

5.494+0.008 

3.1+0.2 

0.82+0.03 

12.6 

E3 

6.289+0.006 

3.1+0.2 

0.74+0.02 

12.5 

E4 

5.638+0.011 

2.2+0.2 

0.71+0.05 

9.2 

El 

6.572+0.006 

3.2+0.3 

0.12+0.03 

8.7 

E2 

- 

- 

- 

- 

E3 

6.295+0.006 

2.8+0.3 

0.90+0.04 

7.6 

E4 

5.685+0.008 

2.3+0.3 

0.87+0.03 

6.2 


Notes. upper half: frequencies found in g’ light curve 
lower half: frequencies found in i’ light curve 

(b) 

errors as given by FAMIAS, more realistic is an error around 0.1-0.2 
mFIz 

reference time: first point of lightcurve 

- g’: BJDtbd = 2456519.5846905 

- i’: BJDtbd = 2456520.53415896 


amplitude in the g’ light curve (F2) was not clearly detected in 
the i’ light curve, but its apparent absence might be explained 
by the high noise level in this data set. The elevated power at 
lower frequencies is likely due to inaccuracies in the binary light 
curve modeling, along with the atmospheric extinction and trans¬ 
parency corrections. A much longer time base is needed to im¬ 
prove the frequency determination and to be able to use the pul¬ 
sations for asteroseismology. Consequently, we do not perform 
a more thorough pulsation analysis than this and limit our result 
simply to the detection of pulsations alone. To prepare the light 
curve for binary modeling, we use the results from Table 1 to 
subtract the detected pulsations from the original light curves. 

4.2. Binary Light Curve Modeiing 

The binary light curve exhibits all the typical features of an 
HW Vir system. Due to the short period and the relatively high 
temperature of the subdwarf, the reflection effect is rather strong. 


with an amplitude around 10%. The secondary eclipse appears 
to be nearly total, in accordance with the high inclination and 
the very deep primary eclipse. As our light curves only cover 
one complete orbital cycle, it is not possible to determine an ac¬ 
curate period from our data alone. For the ephemeris we hence 
cite the period derived by Van Noordet al. (2013), which was 
determined from a much larger baseline. We were able to deter¬ 
mine precise eclipse timings using our data by fitting parabolas 
to the minima. They are summarized in Table 2. Using these val¬ 
ues and the period from Van Noord et al. (2013), the ephemeris 
of the primary minimum is given by 

BJDtbd = 2456519.64027(1) H- 0.065817833(83) ■ E (1) 

where E is the cycle number. The error in the period quoted by 
Van Noord et al. (2013) is likely to be an underestimate given 
the omission of systematics and the poor sampling. We believe 
an error of 0.0001 to be more appropriate. A comparison of 
the secondary eclipse compared to the primary eclipse of both 
lightcurves separately reveals a slight departure of the secondary 
mid-eclipse from phase 0.5; such an offset can be caused by both 
the Rpmer delay (extra light travel time due to the binary orbit) 
and an eccentricity e > 0. For small eccentricities, the total shift 
of the secondary eclipse with respect to phase 0.5 is given in 
Barlow et al. (2012): 

A, A, A, ^^sdB/1 2P 

AtsE — AfR0n;ie]; -H Atgcc “ I ^ I ^ COS CiJ 

7TC \q I 7T 

From both light curves we measure a shift of 3 + 1 s between the 
time of the secondary minimum and phase 0.5. With our system 
parameters we would expect a theoretical shift of the secondary 
eclipse with respect to phase 0.5 due to the Rpmer delay of 2 s. 
If we take that into account, we would get a maximal eccentric¬ 
ity of e ■ cos a> < 0.00055. Ignoring the Rpmer delay results in a 
maximum eccentricity of e ■ cos to < 0.0011. Because of the large 
errors in the radial velocity determinations, the eccentricity can¬ 
not be constrained by the radial velocity curve to that precision. 

A photometric solution to the binary light curve (with pul¬ 
sations removed; Fig. 6) was determined using MORO (Mod¬ 
ified ROche program, see Drechsel et al. 1995). This program 
calculates synthetic light curves which we fit to the observations 
using the SIMPLEX algorithm. This light curve solution code 
is based on the Wilson-Devinney approach (Wilson & Devinney 
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Table 2. Eclipse times 


filter 

eclipse 

BJDTBD[d] 

g’ 

primary 

2456519.64027(1) 

g’ 

secondary 

2456519.60733(2) 

i’ 

primary 

2456520.56176(1) 

i’ 

secondary 

2456520.59463(1) 

i’ 

primary 

2456520.62757(1) 



0.4 0.6 0.8 1 1.2 1.4 1.6 

orbital phase 


Fig. 6. g’ and i’ light curves after the removal of the pulsation signal as 
explained in Sect. 4.1 together with the best-htting light curve model. 
The residuals displayed in the two lower panels still show signs of low- 
amplitude pulsations. 


1971) but uses a modified Roche model that considers the mu¬ 
tual irradiation of hot components in close binary systems. More 
details of the analysis method are described in Schaffenroth et al. 
(2013). To calculate the synthetic light curves, 12 -i-5n (n is the 
number of light curves) parameters are used. Such a high num¬ 
ber of partially-correlated parameters will inevitably cause se¬ 
vere problems if too many and wrong combinations are adjusted 
simultaneously. In particular, there is a strong degeneracy with 
respect to the mass ratio. After the orbital inclination, this param¬ 
eter has the strongest effect on the light curve, and it is highly 
correlated with the component radii. Hence we kept the mass 
ratio fixed at certain values and calculated solutions for these 
mass ratios, which were subsequently compared and evaluated 
according to criteria explained below. 

Given the large number of parameters present in the code, it 
is imperative to constrain as many parameters as possible based 
on independent inputs from spectroscopic analyses or theoreti¬ 
cal constraints. From the spectroscopic analysis, we derived the 
effective temperature and the surface gravity of the sdB primary 
and fixed these parameters during the fitting. Due to the early 
spectral type of the primary star, the gravity darkening exponent 
was fixed at = 1, as expected for radiative outer envelopes 
(vonZeipel 1924). For the cool convective companion, g 2 was 
set to 0.32 (Lucy 1967). The linear limb darkening coefficients 
were extrapolated from the table of Claret & Bloemen (201 1). 


To determine the quality of the light curve fit, the sum of 
squared residuals cr of all observational points with respect to 
the synthetic curve was calculated as a measure of the goodness 
of the fit. Unfortunately, the cr values of the best light curve fits 
for the different mass-ratios did not differ significantly. There¬ 
fore, we cannot determine a unique solution from the light curve 
analysis alone. The full set of parameters describing the best- 
fitting solution for a mass ratio of q' = 0.146, which corresponds 
to an sdB with the canonical mass of 0.47 Mq are given in Table 
3, with errors determined by the bootstrapping method. The light 
curves in the g' and i' bands are displayed in Fig. 6 together with 
the best-fit models for these parameters. The parameters of the 
system derived from this lightcurve solution together with the 
semi-amplitude of the RV curve are summarized in Table 4. The 
errors result from error propagation of the errors in K and P. 


5. The brown dwarf nature of the companion 

From the semi-amplitude of the radial velocity curve, the orbital 
period, and the inclination, we can derive the masses and radii of 
both components for each mass ratio. The masses follow from: 


Ml = 


PK^^ 

2 ^ 


iq + 1 )^ 

(q ■ sin 0^ 


M 2 — q ' M\ 


and the fractional radii of the light curve solution together with 


P Ki 

2n sin i 



yield the radii. For each mass ratio we get different masses and 
radii. It was stated already in Sect. 4.2 that it is not possible to de¬ 
termine the mass ratio from the light curve analysis alone. How¬ 
ever, from the sdB mass and radius determined by the light curve 
analysis we can calculate a photometric surface gravity and com¬ 
pare that to the surface gravity derived by the spectroscopic anal¬ 
ysis. The result is shown in Fig. 7. An agreement of spectro¬ 
scopic and photometric surface gravity values is reached for so¬ 
lutions that result in sdB masses between 0.35 and 0.62 Mq. It is 
therefore possible to find a self-consistent solution. This is a for¬ 
tunate situation, because gravity derived from photometry was 
found to be inconsistent with the spectroscopic result in other 
cases, such as AA Dor (Vuckovic et al. 2008). 

To constrain the solutions even more, we can also use the¬ 
oretical mass-radius relations for the low-mass companions by 
Baraffe et al. (2003) and compare them to the masses and radii 
of the companion derived by the light curve solutions for the 
various mass ratios. This was done in a similar way as in 
Schaffenroth et al. (2014b). This comparison is displayed in Fig. 
8. Relations for different ages of 1, 5 and 10 Gyrs were used. 
The measured mass-radius relation is well matched by theoreti¬ 
cal predictions for stars > 3 Gyrs for companion masses between 
0.056 Mq and 0.073 Mq. The corresponding mass range for the 
sdB star extends from 0.35 Mq to 0.53 Mq. 

However, inflation effects have been found in the case of 
hot Jupiter exoplanets (e.g. Udalski et al. 2008) and also in the 
MS -I- BD binary CoRoT-15b (Bouchy et al. 201 1). As the com¬ 
panion is exposed to intense radiation of a luminous hot star at 
a distance of only 0.56 R©, this effect cannot be neglected and 
would result in an underestimation of the radius, if compared to 
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Table 3. Adopted light curve solution. 


Fixed parameters: 


q{=M2lMi) 


0.146 

Teff(l) 

[K] 

33000 

si 


1.0 

si 


0.32 

xiig'y 


0.21 

XldJ 


0.14 

^2 


0.0 

Adjusted parameters: 


i 

[°] 

86.83 + 0.45 

Teff(2) 

[K] 

2960 + 550 



1.0 + 0.002 

^2 


1.2 + 0.05 



4.10 + 0.05 



2.389 + 0.008 


0.99995 + 0.00007 

(i'Y 


0.99926 + 0.00068 



0.026 + 0.01 

X2(g') 


0.44 + 0.06 

xid') 


0.62 + 0.07 

hig'Y 


0.007 + 0.001 

hd'Y 


0.0 + 0.0 

Roche radii*: 

ri(pole) 

[a] 

0.246 + 0.001 

ri (point) 

[a] 

0.249 + 0.002 

ri(side) 

[a] 

0.249 + 0.002 

ri(back) 

[a] 

0.249 + 0.002 

r2(pole) 

[a] 

0.150 + 0.001 

r2(point) 

[a] 

0.154 + 0.002 

r2(side) 

[a] 

0.152 + 0.001 

r2(back) 

[a] 

0.158 + 0.001 


Notes. 

“ Bolometric albedo 
Gravitational darkening exponent 

^ Linear limb darkening coefficient; taken from Claret & Bloemen 
( 2011 ) 

Radiation pressure parameter, see Drechsel et al. (1995) 

Fraction of third light at maximum 
) Roche potentials 

® Relative luminosity; Li is not independently adjusted, but recomputed 
from r 2 and 7’eff(2) 

* Fractional Roche radii in units of separation of mass centres 


non-irradiated models (Baraffe et al. 2003). The maximum infla¬ 
tion effect can be estimated from the comparison of our solutions 
to the theoretical mass-radius relations shown in Fig. 8. This fig¬ 
ure shows that an inflation of more than about 10% can be ex¬ 
cluded, because otherwise none of the theoretical mass-radius 
relations would match the measured one, even if the star were as 
old as 10 Gyrs. 

The mass-radius relation for the companion would be in per¬ 
fect agreement with the light curve solution for a companion 
with a mass of 0.069 Mq, a radius of O.O86.R0, and an age of 
~ 5-10 Gyrs, if we assume an inflation of 6-11 %. The corre¬ 
sponding mass of the sdB is 0.47 Mq, exactly the canonical sdB 
mass, which we therefore adopt for the sdB throughout the rest 


Table 4. Parameters of V2008-1753 



V2008 

-1753 

coordinates 

20 08 16.355 -17 53 10.52 (J2000.0) 

cv 

[mag] 

16.8 

i 

0 

86.83 + 0.45 

K 

[kms-i] 

54.6 + 2.4 

P 

[h] 

1.5796280 + 0.0000002 

ALsdB 

[Mq\ 

0.47 + 0.03 

-^comp 

[Mq\ 

0.069 + 0.005 

a 

[R<3\ 

0.56 + 0.02 

^sdB 

[R<3\ 

0.138 + 0.006 

^comp 

[R<3\ 

0.086 + 0.004 

log g(sdB, phot) 


5.83+0.02 

log g(sdB, spec) 


5.83+0.05 

7’efr,.sdB 

[K] 

32800 + 750 



0 0.2 0.4 0.6 0.8 1 1.2 1.4 


^sdB [^©l 

Fig. 7. Comparison of the photometric and spectroscopic surface gravity 
for the solutions with different mass ratio q = 0.10, 0.11, 0.12, 0.13, 
0.14, 0.145, 0.146, 0.15, 0.16, 0.17, 0.18, 0.19 (marked by the error 
cross). The spectroscopic surface gravity with uncertainty is given by 
the shaded area. 

of the paper. A similar result was found in the analysis of the 
sdB +BD binary J1622 (Schaffenroth et al. 2014b), which has a 
comparable period and parameters. 

The companion has a mass below the limit for hydrogen¬ 
burning and thus appear to be a brown dwarf - the third con¬ 
firmed around a hot subdwarf star. 


6. Summary and conclusions 

We performed an analysis of the spectrum and light curve of 
V2008-1753 and find that this eclipsing binary consists of a pul¬ 
sating sdB with a brown dwarf companion. This is the first sys¬ 
tem of this kind ever found. Similar to J1622, an inflation of the 
brown dwarf by more than about 10 % can be excluded. 

V2008-1753 has the shortest period of all known HW Vir 
systems and the second shortest period of any sdB binary dis¬ 
covered to date. Due to the small separation distance and high 
temperature of the sdB, the amplitude of the reflection effect 
is relatively large (more than 10%). Consequently, this system 
might provide the chance to detect and evaluate spectral features 
of the irradiated companion, similar to AA Doradus and HW Vir 
(Vuckovic et al. 2008, 2014). If the companion’s spectral fea¬ 
tures are detected, the radial velocities of both components could 
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Fig. 8. Comparison of theoretical mass-radius relations of brown dwarfs 
by Baraffe et al. (2003) for an age of 1 Gyr (filled triangles), 5 Gyrs 
(triangles) and 10 Gyrs (diamond) to results from the lightcure analysis. 
Each error cross represents a solution from the light curve analysis for a 
different mass ratio (q = 0.10, 0.11, 0.12, 0.13, 0.14, 0.145, 0.146, 0.15, 
0.16, 0.17,0.18,0.19). The dashed vertical lines mark different values of 
the corresponding sdB masses. The solid line marks the most probable 
solution with q = 0.146, which results in an sdB mass of 0.47 M^. 


be determined. We then could derive an unbiased mass ratio of 
the system and obtain a unique spectroscopic and photometric 
solution. 

Higher quality and resolution spectra should also allow 
the detection of the Rossiter-McLaughlin (RM) effect (cp. 
Vuckovic et al. 2007b). This effect is due to the selective block¬ 
ing of the light of the rotating star during an eclipse. The ampli¬ 
tude is mainly depending on radius ratio, the rotational velocity 
of the primary star, and the inclination of the system. As our sys¬ 
tem has a high inclination and a high radius ratio, the expected 
amplitude is quite high. From the RM effect it is, hence, possi¬ 
ble to determine the rotational velocity independent from spec¬ 
tral line modelling. With the help of the rotational velocity we 
can check, if the sdB is synchronized, which should be the case 
due to synchronization theories, but is currently under debate 
(Schaffenroth et al. 2014b). 

The presence of stellar pulsations in the sdB offers fur¬ 
ther possibilities to characterize this system. Once the pulsa¬ 
tion frequencies are fully resolved and their modes identified 
with higher-precision data over a longer time base, the mass and 
other properties of the sdB can be constrained by asteroseismoT 
ogy. These results can be compared to the light curve models. 
To date, only NY Vir (Van Grootel et al. 2013) has offered this 
opportunity. However, no signature of the companion could be 
identified in the spectrum of this system, and the masses from 
the light curve analysis remain biased. 

The combined presence of pulsations and eclipses, 
furthermore, offers the possibility of eclipse mapping, 
as done for example for NY Vir (Reed et al. 2005; 
Reed & Whole Earth Telescope Xcov 21 and 23 Collaborations 
2006). Thereby, the eclipses can be used to determine pulsation 
modes, which are often difficult to uniquely identify. During 
the eclipse parts of the star are covered. Changing the amount 
and portion of regions of the star visible affects the pulsation 
amplitudes. The effect is changing for different pulsation modes, 
which can be identified in this way. 

V2008-1753 has the potential to eventually replace NYVir 
as the benchmark system for understanding sdB stars and their 


binary nature, if emission lines of the companion are detected. It 
would permit the direct comparison of independent techniques 
(namely light curve modelling, asteroseismology, spectroscopy, 
and radial velocity variations) used to derive the stellar param¬ 
eters. Hence, the reliability of these methods and models could 
be checked as well as possible systematic errors of the derived 
parameters could be further investigated. Most notably, the mass 
determined by asteroseismology could be checked at a high pre¬ 
cision level, if the semi-amplitudes of the radial velocity curves 
could be determined for both components. 
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